Exploring the Feasibility of 16S rRNA Short Amplicon Sequencing-Based Microbiota Analysis for Microbiological Safety Assessment of Raw Oyster.

16S rRNA short amplicon sequencing-based microbiota profiling has been thought of and suggested as a feasible method to assess food safety. However, even if a comprehensive microbial information can be obtained by microbiota profiling, it would not be necessarily sufficient for all circumstances. To prove this, the feasibility of the most widely used V3-V4 amplicon sequencing method for food safety assessment was examined here. We designed a pathogen (Vibrio parahaemolyticus) contamination and/or V. parahaemolyticus-specific phage treatment model of raw oysters under improper storage temperature and monitored their microbial structure changes. The samples stored at refrigerator temperature (negative control, NC) and those that were stored at room temperature without any treatment (no treatment, NT) were included as control groups. The profiling results revealed that no statistical difference exists between the NT group and the pathogen spiked- and/or phage treated-groups even when the bacterial composition was compared at the possible lowest-rank taxa, family/genus level. In the beta-diversity analysis, all the samples except the NC group formed one distinct cluster. Notably, the samples with pathogen and/or phage addition did not form each cluster even though the enumerated number of V. parahaemolyticus in those samples were extremely different. These discrepant results indicate that the feasibility of 16S rRNA short amplicon sequencing should not be overgeneralized in microbiological safety assessment of food samples, such as raw oyster.

The first reporting of prevalence Vibrio species and expression of HSP genes in rayed pearl oyster (Pinctada radiata) under thermal conditions.

The main goal of the present study was to evaluate the influence of thermal exposure on Vibrio population and HSP genes expression (HSP 90, HSP70, and HSP20) in rayed pearl oyster (P. radiata). To this end, the oysters were reared for 30 days at temperatures of 22 °C (control), 25 °C, 27 °C, and 29 °C. The results showed that five dominate Vibrio strains including Vibrio hepatarius, V. harveyi, V. alginolyticus, V. parahaemolyticus, and V. rotiferianus were identified. The highest population of V. parahaemolyticus, V. alginolyticus, and V. harveyi, was found in 29οC group. According to real-time PCR, mantle exhibited the highest expression levels of HSP20, HSP70, and HSP90 genes. A higher level of HSP20 expression was observed at high temperatures (25 °C, 27 °C, and 29 °C) in the gonad and mantle compared to the control group (22 °C) while decrease in HSP90 expression level was recorded in 25 °C, 27 °C, and 29 °C groups. HSP20 expression level in adductor muscle was remarkably down-regulated in 27 °C and 29 °C groups. In this tissue, HSP70 was detected at highest levels in the 29οC group. In mantle, HSP90 gene expression was lowest at 22 °C water temperature. Several Vibrio strains have been identified from pearl Gulf oyster that haven't been previously reported. The identification of dominant Vibrio species is essential for epidemiological management strategies to control and prevent Vibrio outbreaks in pearl oyster farms. The expression pattern of HSP genes differs in rayed pearl oyster tissues due to differences in their thermal tolerance capability and physiological and biological characteristics. The present study provides useful molecular information for the ecological adaptation of rayed pearl oysters after exposure to different temperature levels.

Early gonadal differentiation is associated with the antagonistic action of Foxl2 and Dmrt1l in the Pacific oyster.

As the second largest phylum in the zoological kingdom next to arthropods, the mechanism of gonadal differentiation in mollusca is quite complex. Currently, although much has been carried out on gonadal differentiation in the Pacific oyster, there is still unknown information that needs to be further explored. Here, analysis of the Foxl2 and Dmrt1l expression in samples at different development periods of male and female gonads as well as in annual gonad samples revealed that Log10 (Foxl2/Dmrt1l) values were an effective method for sex identification in oysters. In differentiated gonadal tissue, Log10 (Foxl2/Dmrt1l) values greater than 2 were females and less than 1 for males. Subsequent sequential sampling of the same individuals verified that Log10 (Foxl2/Dmrt1l) values greater than 2 for resting gonads would develop as females and less than 1 would develop as males in the future. Relative expression analysis of Foxl2 and Dmrt1l in the annual samples revealed a negative correlation between Log10 (Foxl2) and Log10 (Dmrt1l). Double fluorescence reporter validation results showed that DMRT1L protein was able to bind the Foxl2 promoter and repress its activity with a weak dosage effect. Antagonism between Dmrt1l and Foxl2 is therefore not restricted to vertebrates, and the competing regulatory networks are of great significance in the maintenance of gonadal sex in oysters after sexual differentiation. This study provides novel ideas and insights into the study of early gonadal differentiation in the adult oyster.

Migration of repetitive DNAs during evolution of the permanent translocation heterozygosity in the oyster plant (Tradescantia section Rhoeo).

Due to translocation heterozygosity for all chromosomes in the cell complement, the oyster plant (Tradescantia spathacea) forms a complete meiotic ring. It also shows Rabl-arrangement at interphase, featured by polar centromere clustering. We demonstrate that the pericentromeric regions of the oyster plant are homogenized in concert by three subtelomeric sequences: 45S rDNA, (TTTAGGG)n motif, and TSrepI repeat. The Rabl-based clustering of pericentromeric regions may have been an excellent device to combine the subtelomere-pericentromere sequence migration (via inversions) with the pericentromere-pericentromere DNA movement (via whole arm translocations) that altogether led to the concerted homogenization of all the pericentromeric domains by the subtelomeric sequences. We also show that the repetitive sequence landscape of interstitial chromosome regions contains many loci consisting of Arabidopsis-type telomeric sequence or of TSrepI repeat, and it is extensively heterozygous. However, the sequence arrangement on some chromosomal arms suggest segmental inversions that are fully or partially homozygous, a fact that could be explained if the inversions started to create linkages already in a bivalent-forming ancestor. Remarkably, the subterminal TSrepI loci reside exclusively on the longer arms that could be due to sharing sequences between similarly-sized chromosomal arms in the interphase nucleus. Altogether, our study spotlights the supergene system of the oyster plant as an excellent model to link complex chromosome rearrangements, evolution of repetitive sequences, and nuclear architecture.

Teneurin and TCAP Phylogeny and Physiology: Molecular Analysis, Immune Activity, and Transcriptomic Analysis of the Stress Response in the Sydney Rock Oyster (Saccostrea glomerata) Hemocytes.

Teneurin C-terminal associated peptide (TCAP) is an ancient bioactive peptide that is highly conserved in metazoans. TCAP administration reduces cellular and behavioral stress in vertebrate and urochordate models. There is little information for invertebrates regarding the existence or function of a TCAP. This study used the Sydney rock oyster (SRO) as a molluscan model to characterize an invertebrate TCAP, from molecular gene analysis to its physiological effects associated with hemocyte phagocytosis. We report a single teneurin gene (and 4 teneurin splice variants), which encodes a precursor with TCAP that shares a vertebrate-like motif, and is similar to that of other molluscan classes (gastropod, cephalopod), arthropods and echinoderms. TCAP was identified in all SRO tissues using western blotting at 1-2 different molecular weights (~22 kDa and ~37kDa), supporting precursor cleavage variation. In SRO hemolymph, TCAP was spatially localized to the cytosol of hemocytes, and with particularly high density immunoreactivity in granules. Based on 'pull-down' assays, the SRO TCAP binds to GAPDH, suggesting that TCAP may protect cells from apoptosis under oxidative stress. Compared to sham injection, the intramuscular administration of TCAP (5 pmol) into oysters modulated their immune system by significantly reducing hemocyte phagocytosis under stress conditions (low salinity and high temperature). TCAP administration also significantly reduced hemocyte reactive oxygen species production at ambient conditions and after 48 h stress, compared to sham injection. Transcriptomic hemocyte analysis of stressed oysters administered with TCAP demonstrated significant changes in expression of genes associated with key metabolic, protective and immune functions. In summary, this study established a role for TCAP in oysters through modulation of physiological and molecular functions associated with energy conservation, stress and cellular defense.

Detection and Genetic Correlation Analysis of Diarrhea Cases and Norovirus in Oysters in Yantai, China.

Background: This study aimed to assess the correlation between Norovirus (NoV), diarrhea, and raw oysters from the eastern coastal areas of Yantai, Shandong, China. Methods: Marine oysters were selected from the three aquatic markets in Laishan district, Yantai City, in March 2019. Meanwhile, 100 fecal samples were collected from patients with diarrhea from the same areas during the same period. Nucleic acids were extracted from these samples and detected by employing reverse transcription polymerase chain reaction (RT-PCR) for NoV GI/GII. The VP1 gene of the coat protein of NoV was amplified by semi-nested RT-PCR and sequenced. Sequence comparison of VP1 was performed with BioEdit software, and the evolutionary tree was constructed with Mega7.0 software. Results: Of the 151 oysters, 42 (27.8%) were positive for NoV. Among them, 32 (21.2%) were GII-positive, 10 (6.6%) were GI-positive, and one GI VP1 sequence was obtained in the oyster samples. Of 100 fecal samples from patients with diarrhea, 38 were GII-positive and 17 were GI-positive. Totally, 19 GII VP 1 sequences and eight GI VP 1 sequences were obtained. Two G1 VP 1 sequences in two fecal samples showed 98.7% nucleotide sequence identity and 99.1% amino acid sequence identity G1 VP 1 acquired in the oyster sample. Conclusions: The results suggest that oysters may be responsible for the spread of NoV in Yantai, Shandong province, China.

A microcosm study of microbial community profiles during sediment remediation using pyrolyzed oyster shells.

The accumulation of organic and inorganic components in sediments leads to a deterioration in the environment and an imbalance in the coastal ecosystem. Currently, capping is the most effective technology for remediating polluted sediment and restoring ecosystems. A microcosm experiment was designed using pyrolyzed oyster shell (POS). These were mixed in with coastal sediment or added as a capping layer. The results showed that POS effectively decreased pollutants, including PO4-P and NH4-N. Metagenomics analysis was performed using 16S rRNA gene sequencing and the most abundant phyla identified in the POS treated and untreated sediments were Proteobacteria, followed by Firmicutes, Bacteroidetes, Chloroflexi, Fusobacteria, Nitrospirae, and Spirochaetes. The relative abundance of Proteobacteria members of the Class Gammaproteobacteria significantly increased, but Deltaproteobacteria gradually decreased throughout the experiment in POS-covered sediment. This suggests that the POS effectively promoted a shift from anaerobic to facultative anaerobic or aerobic microbial communities in the sediment. Dominant species of facultative anaerobic or microaerophilic bacteria from the order Chromatiales and phylum Nitrospirae were observed in the POS-covered sediment. Based on these study results, it can be concluded that POS is an effective covering material for sediment remediation and restores the microbial communities in sediments.

Variation in Survival and Gut Microbiome Composition of Hatchery-Grown Native Oysters at Various Locations within the Puget Sound.

The Olympia oyster (Ostrea lurida) of the Puget Sound suffered a dramatic population crash, but restoration efforts hope to revive this native species. One overlooked variable in the process of assessing ecosystem health is association of bacteria with marine organisms and the environments they occupy. Oyster microbiomes are known to differ significantly between species, tissue type, and the habitat in which they are found. The goals of this study were to determine the impact of field site and habitat on the oyster microbiome and to identify core oyster-associated bacteria in the Puget Sound. Olympia oysters from one parental family were deployed at four sites in the Puget Sound both inside and outside of eelgrass (Zostera marina) beds. Using 16S rRNA gene amplicon sequencing of the oyster gut, shell, and surrounding seawater and sediment, we demonstrate that gut-associated bacteria are distinct from the surrounding environment and vary by field site. Furthermore, regional differences in the gut microbiota are associated with the survival rates of oysters at each site after 2 months of field exposure. However, habitat type had no influence on microbiome diversity. Further work is needed to identify the specific bacterial dynamics that are associated with oyster physiology and survival rates. IMPORTANCE This is the first exploration of the microbial colonizers of the Olympia oyster, a native oyster species to the West Coast, which is a focus of restoration efforts. The patterns of differential microbial colonization by location reveal microscale characteristics of potential restoration sites which are not typically considered. These microbial dynamics can provide a more holistic perspective on the factors that may influence oyster performance.

Diversification of the aquaporin family in geographical isolated oyster species promote the adaptability to dynamic environments.

BACKGROUND: The diversified aquaporin (AQP) family that was derived from gene duplication and subsequent functional differentiation play critical roles in multiple physiological processes and in adaptation to the dynamic environments during the evolutionary process. Oysters are a group of bivalve fauna in Mollusca that were widely distributed around the world and show extraordinary adaptation to harsh environments. However, knowledge is lacking with the diversity and evolution of the AQP family in oysters, even in molluscs. RESULTS: Here, we performed a comprehensive analysis of the AQP family in three geographical isolated oyster species that are native to different environments. Genome distribution and phylogenetic analysis revealed that the expansion of the AQP family in oysters were attributed to tandem duplication. Synteny analysis indicated that large-scale inversions lead to the independent duplication or deletion of the AQPs after speciation. As a consequence, these independent duplication events contributed to the diversification of the AQP family in different oysters. Pore pattern analysis suggested that the duplicated AQPs in oysters were highly diversified in inner surface profiles, implying the subsequent functional differentiation. The comparison conducted based on the transcriptome data demonstrated that the functional differentiated AQP family members in oysters may play critical roles in maintaining the balance between the stationary homeostasis and dynamic environments. CONCLUSIONS: Our observation provides evidence for the correlation between the duplicated and functional differentiated AQP family and the adaptation to stationary life under dynamic environments in oysters. Additionally, it also broadens our knowledge of the evolution of AQP family in molluscs.

Revealing the composition of the eukaryotic microbiome of oyster spat by CRISPR-Cas Selective Amplicon Sequencing (CCSAS).

BACKGROUND: The microbiome affects the health of plants and animals, including humans, and has many biological, ecological, and evolutionary consequences. Microbiome studies typically rely on sequencing ribosomal 16S RNA gene fragments, which serve as taxonomic markers for prokaryotic communities; however, for eukaryotic microbes this approach is compromised, because 18S rRNA gene sequences from microbial eukaryotes are swamped by contaminating host rRNA gene sequences. RESULTS: To overcome this problem, we developed CRISPR-Cas Selective Amplicon Sequencing (CCSAS), a high-resolution and efficient approach for characterizing eukaryotic microbiomes. CCSAS uses taxon-specific single-guide RNA (sgRNA) to direct Cas9 to cut 18S rRNA gene sequences of the host, while leaving protistan and fungal sequences intact. We validated the specificity of the sgRNA on ten model organisms and an artificially constructed (mock) community of nine protistan and fungal pathogens. The results showed that > 96.5% of host rRNA gene amplicons were cleaved, while 18S rRNA gene sequences from protists and fungi were unaffected. When used to assess the eukaryotic microbiome of oyster spat from a hatchery, CCSAS revealed a diverse community of eukaryotic microbes, typically with much less contamination from oyster 18S rRNA gene sequences than other methods using non-metazoan or blocking primers. However, each method revealed taxonomic groups that were not detected using the other methods, showing that a single approach is unlikely to uncover the entire eukaryotic microbiome in complex communities. To facilitate the application of CCSAS, we designed taxon-specific sgRNA for ~16,000 metazoan and plant taxa, making CCSAS widely available for characterizing eukaryotic microbiomes that have largely been neglected. CONCLUSION: CCSAS provides a high-through-put and cost-effective approach for resolving the eukaryotic microbiome of metazoa and plants with minimal contamination from host 18S rRNA gene sequences. Video Abstract.

Elevated seasonal temperature disrupts prooxidant-antioxidant homeostasis and promotes cellular apoptosis in the American oyster, Crassostrea virginica, in the Gulf of Mexico: a field study.

One of the major impacts of climate change has been the marked rise in global temperature. Recently, we demonstrated that high temperatures (1-week exposure) disrupt prooxidant-antioxidant homeostasis and promote cellular apoptosis in the American oyster. In this study, we evaluated the effects of seasonal sea surface temperature (SST) on tissue morphology, extrapallial fluid (EPF) conditions, heat shock protein-70 (HSP70), dinitrophenyl protein (DNP, an indicator of reactive oxygen species, ROS), 3-nitrotyrosine protein (NTP, an indicator of RNS), catalase (CAT), superoxide dismutase (SOD) protein expressions, and cellular apoptosis in gills and digestive glands of oysters collected on the southern Texas coast during the winter (15 °C), spring (24 °C), summer (30 °C), and fall (27 °C). Histological observations of both tissues showed a notable increase in mucus production and an enlargement of the digestive gland lumen with seasonal temperature rise, whereas biochemical analyses exhibited a significant decrease in EPF pH and protein concentration. Immunohistochemical analyses showed higher expression of HSP70 along with the expression of DNP and NTP in oyster tissues during summer. Intriguingly, CAT and SOD protein expressions exhibited significant upregulation with rising seasonal temperatures (15 to 27 °C), which decreased significantly in summer (30 °C), leaving oysters vulnerable to oxidative and nitrative damage. qRT-PCR analysis revealed a significant increase in HSP70 mRNA levels in oyster tissues during the warmer seasons. In situ TUNNEL assay showed a significant increase in apoptotic cells in seasons with high temperature. These results suggest that elevated SST induces oxidative/nitrative stress through the overproduction of ROS/RNS and disrupts the antioxidant system which promotes cellular apoptosis in oysters.

Identification of Distant Regulatory Elements Using Expression Quantitative Trait Loci Mapping for Heat-Responsive Genes in Oysters.

Many marine ectotherms, especially those inhabiting highly variable intertidal zones, develop high phenotypic plasticity in response to rapid climate change by modulating gene expression levels. Herein, we examined the regulatory architecture of heat-responsive gene expression plasticity in oysters using expression quantitative trait loci (eQTL) analysis. Using a backcross family of Crassostrea gigas and its sister species Crassostrea angulata under acute stress, 56 distant regulatory regions accounting for 6-26.6% of the gene expression variation were identified for 19 heat-responsive genes. In total, 831 genes and 164 single nucleotide polymorphisms (SNPs) that could potentially regulate expression of the target genes were screened in the eQTL region. The association between three SNPs and the corresponding target genes was verified in an independent family. Specifically, Marker13973 was identified for heat shock protein (HSP) family A member 9 (HspA9). Ribosomal protein L10a (RPL10A) was detected approximately 2 kb downstream of the distant regulatory SNP. Further, Marker14346-48 and Marker14346-85 were in complete linkage disequilibrium and identified for autophagy-related gene 7 (ATG7). Nuclear respiratory factor 1 (NRF1) was detected approximately 3 kb upstream of the two SNPs. These results suggested regulatory relationships between RPL10A and HSPA9 and between NRF1 and ATG7. Our findings indicate that distant regulatory mutations play an important role in the regulation of gene expression plasticity by altering upstream regulatory factors in response to heat stress. The identified eQTLs provide candidate biomarkers for predicting the persistence of oysters under future climate change scenarios.

Reconstruction of the evolutionary biogeography reveal the origins and diversification of oysters (Bivalvia: Ostreidae).

Oysters (Bivalvia: Ostreidae Rafinesque, 1815) live in the intertidal and shallow subtidal areas worldwide. Despite their long evolutionary histories, abundant fossil records, global distribution, and ecological significance, a systematic time-dependent biogeographical analysis of this family is still lacking. Using combined mitochondrial (COI and 16S rRNA) and nuclear (18S rRNA, 28S rRNA, H3 and ITS2) gene makers for 80% (70/88) of the recognized extant Ostreidae, we reconstructed the global phylogenetic and biogeographical relationships throughout the evolutionary history of oysters. The result provided a holistic view of the origin, migration and dispersal patterns of Ostreidae. The phylogenetic results and fossil evidence indicated that Ostreidae originated from the circum-Arctic region in the Early Jurassic. The widening of the Atlantic Ocean and changes in the Tethys Ocean further facilitated their subsequent diversification during the Cretaceous and the Palaeogene periods. In particular, Crassostrea and Saccostrea exhibited relatively low dispersal abilities and their major diversifications were consistent with the tectonic events. Environmental adaptations and reproductive patterns, therefore, should play key roles in the formation of oyster distribution patterners, rather than the dispersal ability of their planktonic larvae. The diversity dynamics inferred by standard phylogenetic are consistent with the fossil record, however, further systematic classification, especially for fossil genus Ostrea, would enhance our understanding on extant and fossil oysters. The present study of the historical biogeography of oysters provides new insights into the evolution and speciation of oysters. Our findings also provide a foundation for the assessment of evolutionary patterns and ecological processes in intertidal and inshore life.

Evidence from oyster suggests an ancient role for Pdx in regulating insulin gene expression in animals.

Hox and ParaHox genes encode transcription factors with similar expression patterns in divergent animals. The Pdx (Xlox) homeobox gene, for example, is expressed in a sharp spatial domain in the endodermal cell layer of the gut in chordates, echinoderms, annelids and molluscs. The significance of comparable gene expression patterns is unclear because it is not known if downstream transcriptional targets are also conserved. Here, we report evidence indicating that a classic transcriptional target of Pdx1 in vertebrates, the insulin gene, is a likely direct target of Pdx in Pacific oyster adults. We show that one insulin-related gene, cgILP, is co-expressed with cgPdx in oyster digestive tissue. Transcriptomic comparison suggests that this tissue plays a similar role to the vertebrate pancreas. Using ATAC-seq and ChIP, we identify an upstream regulatory element of the cgILP gene which shows binding interaction with cgPdx protein in oyster hepatopancreas and demonstrate, using a cell culture assay, that the oyster Pdx can act as a transcriptional activator through this site, possibly in synergy with NeuroD. These data argue that a classic homeodomain-target gene interaction dates back to the origin of Bilateria.

Prediction Accuracies of Genomic Selection for Nine Commercially Important Traits in the Portuguese Oyster (Crassostrea angulata) Using DArT-Seq Technology.

Genomic selection has been widely used in terrestrial animals but has had limited application in aquaculture due to relatively high genotyping costs. Genomic information has an important role in improving the prediction accuracy of breeding values, especially for traits that are difficult or expensive to measure. The purposes of this study were to (i) further evaluate the use of genomic information to improve prediction accuracies of breeding values from, (ii) compare different prediction methods (BayesA, BayesCπ and GBLUP) on prediction accuracies in our field data, and (iii) investigate the effects of different SNP marker densities on prediction accuracies of traits in the Portuguese oyster (Crassostrea angulata). The traits studied are all of economic importance and included morphometric traits (shell length, shell width, shell depth, shell weight), edibility traits (tenderness, taste, moisture content), and disease traits (Polydora sp. and Marteilioides chungmuensis). A total of 18,849 single nucleotide polymorphisms were obtained from genotyping by sequencing and used to estimate genetic parameters (heritability and genetic correlation) and the prediction accuracy of genomic selection for these traits. Multi-locus mixed model analysis indicated high estimates of heritability for edibility traits; 0.44 for moisture content, 0.59 for taste, and 0.72 for tenderness. The morphometric traits, shell length, shell width, shell depth and shell weight had estimated genomic heritabilities ranging from 0.28 to 0.55. The genomic heritabilities were relatively low for the disease related traits: Polydora sp. prevalence (0.11) and M. chungmuensis (0.10). Genomic correlations between whole weight and other morphometric traits were from moderate to high and positive (0.58-0.90). However, unfavourably positive genomic correlations were observed between whole weight and the disease traits (0.35-0.37). The genomic best linear unbiased prediction method (GBLUP) showed slightly higher accuracy for the traits studied (0.240-0.794) compared with both BayesA and BayesCπ methods but these differences were not significant. In addition, there is a large potential for using low-density SNP markers for genomic selection in this population at a number of 3000 SNPs. Therefore, there is the prospect to improve morphometric, edibility and disease related traits using genomic information in this species.

Climate change alters the haemolymph microbiome of oysters.

The wellbeing of marine organisms is connected to their microbiome. Oysters are a vital food source and provide ecological services, yet little is known about how climate change such as ocean acidification and warming will affect their microbiome. We exposed the Sydney rock oyster, Saccostrea glomerata, to orthogonal combinations of temperature (24, 28 °C) and pCO2 (400 and 1000 µatm) for eight weeks and used amplicon sequencing of the 16S rRNA (V3-V4) gene to characterise the bacterial community in haemolymph. Overall, elevated pCO2 and temperature interacted to alter the microbiome of oysters, with a clear partitioning of treatments in CAP ordinations. Elevated pCO2 was the strongest driver of species diversity and richness and elevated temperature also increased species richness. Climate change, both ocean acidification and warming, will alter the microbiome of S. glomerata which may increase the susceptibility of oysters to disease.

Reconstruction of ancient homeobox gene linkages inferred from a new high-quality assembly of the Hong Kong oyster (Magallana hongkongensis) genome.

BACKGROUND: Homeobox-containing genes encode crucial transcription factors involved in animal, plant and fungal development, and changes to homeobox genes have been linked to the evolution of novel body plans and morphologies. In animals, some homeobox genes are clustered together in the genome, either as remnants from ancestral genomic arrangements, or due to coordinated gene regulation. Consequently, analyses of homeobox gene organization across animal phylogeny provide important insights into the evolution of genome organization and developmental gene control, and their interaction. However, homeobox gene organization remains to be fully elucidated in several key animal ancestors, including those of molluscs, lophotrochozoans and bilaterians. RESULTS: Here, we present a high-quality chromosome-level genome assembly of the Hong Kong oyster, Magallana hongkongensis (2n = 20), for which 93.2% of the genomic sequences are contained on 10 pseudomolecules (~ 758 Mb, scaffold N50 = 72.3 Mb). Our genome assembly was scaffolded using Hi-C reads, facilitating a larger scaffold size compared to the recently published M. hongkongensis genome of Peng et al. (Mol Ecol Resources, 2020), which was scaffolded using the Crassostrea gigas assembly. A total of 46,963 predicted gene models (45,308 protein coding genes) were incorporated in our genome, and genome completeness estimated by BUSCO was 94.6%. Homeobox gene linkages were analysed in detail relative to available data for other mollusc lineages. CONCLUSIONS: The analyses performed in this study and the accompanying genome sequence provide important genetic resources for this economically and culturally valuable oyster species, and offer a platform to improve understanding of animal biology and evolution more generally. Transposable element content is comparable to that found in other mollusc species, contrary to the conclusion of another recent analysis. Also, our chromosome-level assembly allows the inference of ancient gene linkages (synteny) for the homeobox-containing genes, even though a number of the homeobox gene clusters, like the Hox/ParaHox clusters, are undergoing dispersal in molluscs such as this oyster.

Tracing key genes associated with the Pinctada margaritifera albino phenotype from juvenile to cultured pearl harvest stages using multiple whole transcriptome sequencing.

BACKGROUND: Albino mutations are commonly observed in the animal kingdom, including in bivalves. In the black-lipped pearl oyster Pinctada margaritifera, albino specimens are characterized by total or partial absence of colouration resulting in typical white shell phenotype expression. The relationship of shell colour with resulting cultured pearl colour is of great economic interest in P. margaritifera, on which a pearl industry is based. Hence, the albino phenotype provides a useful way to examine the molecular mechanisms underlying pigmentation. RESULTS: Whole transcriptome RNA-sequencing analysis comparing albino and black wild-type phenotypes at three stages over the culture cycle of P. margaritifera revealed a total of 1606, 798 and 187 differentially expressed genes in whole juvenile, adult mantle and pearl sac tissue, respectively. These genes were found to be involved in five main molecular pathways, tightly linked to known pigmentation pathways: melanogenesis, calcium signalling pathway, Notch signalling pathway, pigment transport and biomineralization. Additionally, significant phenotype-associated SNPs were selected (N = 159), including two located in the Pif biomineralization gene, which codes for nacre formation. Interestingly, significantly different transcript splicing was detected between juvenile (N = 1366) and adult mantle tissue (N = 313) in, e.g., the tyrosinase Tyr-1 gene, which showed more complex regulation in mantle, and the Notch1 encoding gene, which was upregulated in albino juveniles. CONCLUSION: This multiple RNA-seq approach provided new knowledge about genes associated with the P. margaritifera albino phenotype, highlighting: 1) new molecular pathways, such as the Notch signalling pathway in pigmentation, 2) associated SNP markers with biomineraliszation gene of interest like Pif for marker-assisted selection and prevention of inbreeding, and 3) alternative gene splicing for melanin biosynthesis implicating tyrosinase.

Sequence Composition Underlying Centromeric and Heterochromatic Genome Compartments of the Pacific Oyster Crassostrea gigas.

Segments of the genome enriched in repetitive sequences still present a challenge and are omitted in genome assemblies. For that reason, the exact composition of DNA sequences underlying the heterochromatic regions and the active centromeres are still unexplored for many organisms. The centromere is a crucial region of eukaryotic chromosomes responsible for the accurate segregation of genetic material. The typical landmark of centromere chromatin is the rapidly-evolving variant of the histone H3, CenH3, while DNA sequences packed in constitutive heterochromatin are associated with H3K9me3-modified histones. In the Pacific oyster Crassostrea gigas we identified its centromere histone variant, Cg-CenH3, that shows stage-specific distribution in gonadal cells. In order to investigate the DNA composition of genomic regions associated with the two specific chromatin types, we employed chromatin immunoprecipitation followed by high-throughput next-generation sequencing of the Cg-CenH3- and H3K9me3-associated sequences. CenH3-associated sequences were assigned to six groups of repetitive elements, while H3K9me3-associated-ones were assigned only to three. Those associated with CenH3 indicate the lack of uniformity in the chromosomal distribution of sequences building the centromeres, being also in the same time dispersed throughout the genome. The heterochromatin of C. gigas exhibited general paucity and limited chromosomal localization as predicted, with H3K9me3-associated sequences being predominantly constituted of DNA transposons.

The shell matrix of the european thorny oyster, Spondylus gaederopus: microstructural and molecular characterization.

Molluscs, the largest marine phylum, display extraordinary shell diversity and sophisticated biomineral architectures. However, mineral-associated biomolecules involved in biomineralization are still poorly characterised. We report the first comprehensive structural and biomolecular study of Spondylus gaederopus, a pectinoid bivalve with a peculiar shell texture. Used since prehistoric times, this is the best-known shell of Europe's cultural heritage. We find that Spondylus microstructure is very poor in mineral-bound organics, which are mostly intercrystalline and concentrated at the interface between structural layers. Using high-resolution liquid chromatography tandem mass spectrometry (LC-MS/MS) we characterized several shell protein fractions, isolated following different bleaching treatments. Several peptides were identified as well as six shell proteins, which display features and domains typically found in biomineralized tissues, including the prevalence of intrinsically disordered regions. It is very likely that these sequences only partially represent the full proteome of Spondylus, considering the lack of genomics data for this genus and the fact that most of the reconstructed peptides do not match with any known shell proteins, representing consequently lineage-specific sequences. This work sheds light onto the shell matrix involved in the biomineralization in spondylids. Our proteomics data suggest that Spondylus has evolved a shell-forming toolkit, distinct from that of other better studied pectinoids - fine-tuned to produce shell structures with high mechanical properties, while limited in organic content. This study therefore represents an important milestone for future studies on biomineralized skeletons and provides the first reference dataset for forthcoming molecular studies of Spondylus archaeological artifacts.